ENUMERATION OF SNAKES AND 
CYCLE-ALTERNATING PERMUTATIONS 



MATTHIEU JOSUAT- VERGES 
To the memory of Vladimir Arnol'd 

Abstract. Springer numbers are an analog of Euler numbers for the group of signed 
<^^> permutations. Arnol'd showed that they count some objects called snakes, that 

generalize alternating permutations. Hoffman established a link between Springer 
' numbers, snakes, and some polynomials related with the successive derivatives of 

Q , trigonometric functions. 

■ The goal of this article is to give further combinatorial properties of derivative 

polynomials, in terms of snakes and other objects: cycle-alternating permutations, 
C^*) , weighted Dyck or Motzkin paths, increasing trees and forests. We obtain the gen- 

erating functions, in terms of trigonometric functions for exponential ones and in 
terms of J-fractions for ordinary ones. We also define natural g-analogs, make a link 
with normal ordering problems and combinatorial theory of differential equations. 
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1. Introduction 
It is well-known that the Euler numbers E n defined by 



oo _ 
Z 



*i> ' E n — = tan z + sec z (1) 

K : 

count alternating permutations in <3 n , i.e. a such that o\ > 02 < 03 > ...cr n . The 
study of these, as well as other classes of permutations counted by E n , is a vast topic 
in enumerative combinatorics, see the survey of Stanley [21]. By a construction due 
to Springer [20], there is an integer K(W) defined for any Coxeter group W such 
that K(6 n ) = E n . As for the groups of signed permutations, Arnol'd [1] introduced 
some particular kind of signed permutations called snakes, counted by the numbers 
K(&%) and K{&%). Thus snakes can be considered as a "signed analog" of alternating 
^ permutations. Algorithmically, Arnol'd [1] gives a method to compute these integers 

with recurrences organized in triangular arrays similar to the Seidel-Entriger triangle 
of Euler numbers E n (see for example [11]). 

In this article, we are mostly interested in the numbers S n = K(&^), which were 
presviously considered by Glaisher [9, §§ 109 and 119] in another context. Springer [20] 
shows that they satisfy 

"a" 



n! cosz — smz 



! 

n=0 

and we call S n the nth Springer number (although there is in theory a Springer number 
associated with each Coxeter group, this name without further specification usually 
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refer to type B). Another link between snakes and trigonometric functions has been 
established by Hoffman [13], who studies polynomials associated with the successive 
derivatives of tan and sec, i.e. P n (t) and Q n (t) such that: 
d n d n 

-— - tan x = P n (tan x) , -— ■ secx = Q n (tan x) secx. (3) 

Among various other results, Hoffman proves that Q n (l) = S n , and that -Pn(l) = 2 n £' ra 
counts some objects called /3-snakes that are a superset of type B snakes. A first 
question we can ask is to find the meaning of the parameter t in snakes. 

The main goal of this article is to give combinatorial models of these derivative 
polynomials P n and Q n (as well as another sequence Qn related the with the ath 
power of sec), in terms of several objects: 

• Snakes (see Definition 3.2). Besides type B snakes of Arnol'd and /3-snakes of 
Hoffman [13], we introduce another variant. 

• Cycle-alternating permutations (see Definition 3.7). These are essentially the image 
of snakes via Foata's fundamental transform [17]. 

• Weighted Dyck prefixes and weighted Motzkin paths (see Section 3). These are 
two different generalizations of some weighted Dyck paths counted by E n . 

• Increasing trees and forests (see Section 4). The number i?2n+i counts increasing 
complete binary trees, and E^ n "almost complete" ones, our trees and forests are 
a generalization of these. 

Various exponential and ordinary generating functions are proved combinatorially us- 
ing these objects. They will present a phenomenon similar to the case of Euler numbers 
E n : they are given in terms of trigonometric functions for the exponential generating 
functions, and in terms of continued fractions for the ordinary ones. For example, 
in Theorem 3.17, we prove bijectively a continued fraction for the (formal) Laplace 
transform of (cosz — isinz)~ a . While this result was known analytically since long 
ago [22], we show here that it can be fully understood on the combinatorial point of 
view. 



This article is organized as follows. In Section 2, we give first properties of the 
derivative polynomials, such as recurrence relations, generating functions, we also 
introduce a natural o-analog and link this with the normal ordering problem. In 
Section 3, we give combinatorial models of the derivative polynomials in terms of 
subset of signed permutations, and obtain the generating functions. In Section 4, we 
give combinatorial models of the derivative polynomials in terms of increasing trees 
and forests in two different ways: via the combinatorial theory of differential equations 
and via the normal ordering problem. 



2. Derivative polynomials 
2.1. Definitions. Generalizing Q n {t) defined in (3), we consider Qn(t) such that 

d n 

sec a x = (tan x) sec a x. (4) 



dx 

Besides Q n which is the case a = 1 in , of particular interest will be the case a = 2 
denoted by R n = Q n . Since sec 2 x is the derivative of tanx, it follows that there is 
the simple relation P n+ \{t) = (1 + t 2 )R n (t). 
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By differentiation of (3) and (4), one obtains the recurrence relations: 

P n+1 = (1 + t 2 )P n , Qi a ]l = (1 + * 2 )^ Q) + at Qn, (5) 

together with Pq = t and = 1. Let us give other elementary properties which are 
partly taken from [13]. 

Proposition 2.1. The polynomial has the same parity as n, i.e. Qn (—t) = 
(—l) n Q^\t), and P n has different parity. We have: 

P n (l) = 2 n E n , Q n (l) = S n , R n (l) = 2 n E n+l . (6) 

(Recall that Q n = Qn o,nd R n = Qn ■) The generating functions are: 

Z^_ = sinZ + tCOSZ £. = 1 ^ 

^ ™ w n! cosz-tsinz' w , (cos z - t sin z) a ' w 

n=0 n=0 v 7 

Proof. The generating function of {Qn }n>o is n °t present in the reference [13], but 
can be obtained in the same way as in the particular case a = 1 (which is in [13]). By 
a Taylor expansion and a trigonometric addition formula, we have: 

EQl a) (tann)sec a u^- = sec a (it + z) = SeC — , (8) 
^ n K ' n! y ' (cosz - tanitsinz) a v ; 

n=0 v ; 

We can divide on both sides by sec a u, let t = tan it, and the result follows. I 

These exponential generating functions will be obtained combinatorially in Section 3 
using classes of signed permutations. It is interesting to note that P(z, t) = ^2 P n z n /n\ 
is a Mobius transformation as a function of t, in such a way that z i— > (t i— > P(z,t)) 
is a group homomorphism from R to the elliptic Mobius transformations of C fixing i 
and —i. This has an explanation through the fact that 

P(z, t) = tan(arctan t + z), (9) 
and consequently, as observed in [13]: 

P(z,P(z',t)) = P(z + z',t), (10) 
which is the concrete way to say that z i— > (t t— >■ P(z,t)) is a group homomorphism. 

2.2. Operators and ^-analogs. Let D and U be operators acting on polynomials in 
the variable t by 

D(t n ) = [n] q t n -\ U(t n ) = t n+l , (11) 
where [n] q = -j^L _ The first one is known as the q-derivative or Jackson derivative, 
and an important relation is DU — qU D = I, where I is the identity. 

Definition 2.2. Our g-analogs of the polynomials Q n and R n are defined by: 

Q n (t,q) = (D + UDU) n l, R n (t,q) = (D + DUU) n l. (12) 

Of course, 1 should be seen as t° since the operators act on polynomials in t. 

Proposition 2.3. We have 

Qn(-t, q) = (-l) n Q n (t, q), Rn(~t, q) = {-l) n Rn{t, q), (13) 

moreover Q n (t, q) and R n (t, q) are polynomials in t and q with nonnegative coefficients 
such that Q n (t, 1) = Q n (t) and R n (t, 1) = R n {t). 
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Proof. From DU — qUD = /, we can write 



D + UDU = (I + qU 2 )D + U, D + DUU = (I + q U )D + (1 + q)U, (14) 



then from (12) and (14) it follows that 

Qn+l{t, q) = {l+ qt 2 )D(Q n (t, q)) + tQ n (t, q), 
Rn+l{t,q) = (1 + q 2 t 2 )D(R n (t,q)) + (1 + q)tR n (t, 



(15) 
(16) 



which generalize the recurrences for Q n and R n (a = 1 and a = 2 in (5)). The 
elementary properties given in the proposition can be proved recursively using (15) 
and (16). I 

The fact that Q n (t, q) and R n (t, q) have simple forms in terms of D and U makes a 
link with normal ordering [2]. Given an expression f{D,U) of D and U, the problem 
is to find some coefficients c$j such that we have the normal form: 



f(D.U) 



E 

ij>0 



(17) 



When f(D, U) is polynomial, using the commutation relation DU — qUD = I we can 
always find such coefficients and only finitely of them are non-zero. For example, this 
is what we have done in (14). In Section 4, we will use some general results on normal 
ordering to give combinatorial models of P n and Q n (when q = 1). 

What is interesting about normal ordering is that having the coefficients Cij give 
more insight on the quantity f(D,U). For example, consider the case of Q n (t,q) and 
let f(D, U) = {D + UDU) n , then we have 



Q n (t,q) = (D + UDU) n l 



E 

i,3>0 



i>0 



(18) 



i.e. we can directly obtain the coefficients of Q n (t,q) through this normal form. 



2.3. Continued fractions. Let D\ and U\ be the matrices of operators D and U in 
the basis {t l }j e N. Moreover, let W\ be the row vector (t 1 )^ and V\ be the column 
vector (5io)ieN where we use Kronecker's delta. We have W\ = (l,t,t 2 , . . . ), and: 



Z>i 



\(0) 



[2], 




(0)\ 



7 



/o 



\(0) 





1 



(0)\ 



7 



Vi 






V7 



(19) 



From the previous definitions, it follows that: 

Q n (t,q) = W 1 (D 1 + U 1 D 1 U 1 ) n V 1 , R n (t,q) = W 1 (D 1 + D 1 U 1 U 1 ) n V 1 , 
DxUi - qU x Di = /, WWi = tW x , DxVx = 0. 



(20) 
(21) 



The point of writing this is that from (20) we only need the relations in (21) to 
calculate Q n (t,q) and R n (t,q). Even more, if D2, U2, W2, and V2 are a second set 
of matrices and vectors satisfying relations similar to (21), we also have Q n (t,q) = 
W 2 (D 2 + U2D2U2TV2 and R n (t,q) = W 2 (D 2 + D2U2U2TV2. Indeed, when we have 
the relation (17) for a given f(D, U), the coefficients Cj are obtained only using the 
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commutation relation, so the same identity also holds with either (Di, Ui) or (D 2 , U 2 ). 
We can take W 2 = (1,0, 0, . . . ), V 2 = Vt, D 2 = £>i, and: 

/ * (0)\ 
1 tq 



lh 



1 tq 2 



(22) 



\(o) •. ■,/ 

Then we have indeed D 2 U 2 — qU 2 D 2 = I, W 2 U 2 = tW 2 , and D 2 V 2 = 0. What is nice 
about this second set of matrices and vectors is that it enables us to make a link with 
continued fractions. 



Definition 2.4. For any two sequences {bh}h>o and {Xh}h>i, let 3(bh,^h) denote 

1 

Ai z 



3{b h , \ h ) = — . (23) 



1 - b z 



1 h ^ 
1 — b\z 



We will always use z as variable and h as index of the two sequences so that there 
should be no ambiguity in the notation. 

These are called J- fractions, or Jacobi continued fractions. They are linked with 
moments of formal orthogonal polynomials, but let us give the results for Q n (t, q) and 
R n (t,q) before more details about this. 

Proposition 2.5. We have: 

00 00 

J2 Qn{t, q)z n = Z{b%, X Q h ), J2 = Af ), (24) 



n=0 n=0 



where: 



K = tq h ([h) q + [h+ 1],), b£ = tq h (l + q)[h + 1]„ (25) 

A? = (1 + t^ h - X )[h\l A* = (1 + tV h )[h] q [h + l] ff . (26) 

Proof. First, observe that for any matrix M = (rriij)^^, the product W 2 M n V 2 is 
the upper-left coefficient (M n )o,o of M n . We can obtain this coefficient W 2 M n V 2 the 
following way: 

W 2 M n V 2 = rao,ii«tii,i 2 • • • m in _ 2) i n _ 1 m in _ lj0 . (27) 

ii,...,i n _i>0 

When the matrix is tridiagonal, we can restrict the sum to indices such that < 
1, so that these indices are the successive heights in a Motzkin path. Then (27) shows 
that W 2 M n V 2 is the generating function of Motzkin paths of length n with some 
weights given by the coefficients of M, and using a classical argument [8] it follows 
t^EZ^(W2M n V 2 )z n = Z(m h)h ,m h - 1 , h m h , h - 1 ). 

In the present case, it suffices to check that D 2 + U 2 D 2 U 2 and D 2 + D 2 U 2 U 2 are 
tridiagonal and calculate explicitly their coefficients to obtain the result. See also (34) 
and (36) below for a more general result. I 

In particular, it follows from the previous proposition that we have: 

Q 2n (0,q) = E 2n (q), R 2n {0,q) = E 2n+1 (q), (28) 



6 



MATTHIEU JOSUAT- VERGES 



where E2 n (q) and 2?2n+i(<z) are respectively the g-secant and g-tangent numbers de- 
fined by Han, Randrianarivony, Zeng [12] using continued fractions. 

Remark 2.6. From the fact that Q n (l) = S n and i? n (l) = 2 n E n+ i, the previous 
proposition implies 

oo oo 

]T S n z n = 3{2h + 1, 2h 2 ), Yl E ^ = ^ h + M ¥ ii )- ( 29 ) 

n=0 n=0 

We found only one reference mentioning the latter continued fraction, Sloane's OEIS 
[19], but there is little doubt it can be proved by classical methods. For example, 
using a theorem of Stieltjes and Rogers [10, Theorem 5.2.10], the continued fraction 
can probably be obtained through an addition formula satisfied by the exponential 
generating function of {-E n+ i} n >o, this function being the derivative of tanz + secz, 
explicitly (1 — sin z) . Combinatorially, the result can be proved using Andre trees [7, 
Section 5] and the bijection of Frangon and Viennot [10, Chapter 5]. More generally, 
this shows that if we add a parameter x, then Z{h + 1, X — - ) is the generating 
function of the Andre polynomials defined by Foata and Schiitzenberger [7]. 

An important property of J-fractions is the link with moments of (formal) orthog- 
onal polynomials, and we refer to [8, 23] for the relevant combinatorial facts. We con- 
sider here the continuous dual q-Hahn polynomials p n (x;a,b,c\q), or p n (x) for short. 
In terms of the Askey- Wilson polynomials which depend on one other parameter d, 
p n {x) is just the specialization d = (see [15] for the definitions of these classical 
sequences, but our notations differ, in particular because of a rescaling x — > x/2). We 
have the three-term recurrence relation 

xp n (x) = p n +i{x) + (a + \ - A n - C n ) p n {x) + ,4 n _iC n p n _i (x) , (30) 

together with p_i(x) = and po(x) = 1, where 

A n = 1(1 - abq n )(l - acq 11 ) and C n = a(l - q n )(l - bcq n - x ). (31) 

The p n (x) are orthogonal with respect to the scalar product (/, g) \— > L(fg) where L 
is the linear form such that 

oo 

L{x n )z n =z{a+\-A h - C h ,A h ^C h ) . (32) 

n=0 

The quantity L(x n ) is called the nth moment of the orthogonal sequence {p n (x)} n >o- 
Let fj> n (a,b,c) denote this nth moment of the continuous dual g-Hahn polynomials 
p n (x). An elementary calculation shows the following: 

Proposition 2.7. We have 

n ,. \ Mh/gr-V^) , v Hn(iq,-iq,t) , , 

Qnfo ^ = (TT^ ' ?) = (1 _ g )n • ( 33 ) 

Proof. We just have to identify the ordinary generating functions of both sides in each 
identity. This is possible because we know in each case the explicit form of the J- 
fraction expansion, from (32) and Proposition 2.5. See also (34) and (35) below for a 
more general result. I 

This very simple link between our polynomials and the moments /i n (o; b, c) is one of 
the properties indicating that Q n {t,q) and R n (t,q) are interesting g-analogs of Q n (t) 
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and R n (t). Note that in view of (15) and (16), it is tempting to define a (/-analog 
Qn^itil) °f Qn"\t) by the recurrence: 

QitlM = (l + q a t 2 )D(Q^(t,q)) + [a] q tQ^(t,q), (34) 

together with Q^\t, q) = 1. Though we will not study these apart the particular cases 
a = 1 and a = 2, it is worth mentioning that a natural generalization of the previous 
proposition holds, more precisely we have: 

Q^) = ^Xip t \ (35) 

which can be proved by calculating explicitly the matrix (/ + q a [/| )Z?2 + [ a ]qU2- The 
entries in this matrix show that the generating function of {Qn(t, q)} has a J-fraction 
expansion with coefficients: 

b h = tq h [[h-l + a] q + q a - l [h + l] q ) , X h = (1 + t 2 q 2h ~ 2+a )[h] q [h - 1 + a] q . (36) 

Indeed we have the same coefficients if we replace (a, b, c) with (iq a / 2 , —iq a ^ 2 ,t) in (31) 
and (32). 

3. Combinatorial models via signed permutations 

Definition 3.1. We will denote [n] = {1 . . . n} and [nj = {— n . . . —1} U {1 . . . n}. A 

signed permutation is a permutation ir of [n] such that 7r(— i) = —ir(i) for any i G [n]. 
It will be denoted ir = tt\ . . . ir n where 7Tj = tt(i), and 6^ is the set of all such ir. We 
will also use the cycle notation, indicated by parenthesis: for example ir = 3, —1,2,4 
is (1, 3, 2, — 1, — 3, — 2)(4)(— 4) in cycle notation. 

We consider here two classes of signed permutations: snakes and cycle-alternating 
permutations. These will give combinatorial models of the derivative polynomials. As 
for the g-analog, going through other objects (weighted Dyck prefixes) we can obtain 
a statistic on cycle-alternating permutations counted by the parameter q. 

3.1. Snakes. There are several types of snakes to be distinguished, so let us first give 
the definition: 

Definition 3.2. A signed permutation ir = tti . . . ir n is a snake if tti > tt2 < tt^ > 
...7r n . We denote by S n C &^ the set of snakes of size n. Let 5° C S n be the 
subset of 7r satisfying tt\ > 0, and S®° C 5° be the subset of it satisfying i\\ > and 
(-l)*7r n <0. 

The objects introduced by Arnol'd [1] are the ones in 5°, and are usually called 
snakes of type B. He shows in particular that the cardinal of is the nth Springer 
number S n . The objects in S n are called /3-snakes by Hoffman [13], and he showed 
that the cardinal of this set is -P n (l) = 2 n E n . 

It is practical to take a convention for the values 7To and 7r n+ i. If it G S n (respectively, 
i G 5jJ, 7T S 5°°) we define ttq = — (n + 1) and 7r n+ i = (— l) n (n + 1) (respectively, 
7To = and ir n+ i = ( — l) n (n + 1), ttq = TT n+ i = 0). This is not consistent with 
the inclusions 5°° C 5° C S n , so we can think of these sets to be disjoint to avoid 
ambiguity in the sequel. Then an element of any of these three sets is a tt G suc h 
that ttq < 7Ti > 7T2 < VT3 > . . .7r n+ i. For example in the case of 5^°, the condition 
(— l) n 7r n < means that if n is odd, 7r n > 7r n+ i = 0, and if n is even, 7r n < ir n+ i = 0. 
For clarity, we will often write 7r = (vtq), n%, . . . , 7r n , (7r n+ i). 
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Definition 3.3. Let tt be a snake in one of the sets S n , S®, or S° , and ttq, 7r n+ i as 
above. We define a statistic sc(7r) as the number of sign changes through the values 
7T . . . 7r n+ i, i.e. sc(7r) = #{ j I < i < n, TTiiTi+i < }. 

Using this statistic, we have combinatorial models of the derivative polynomials. 

Theorem 3.4. For any n > 0, we have: 

Pn(t) = tSC(7r) < Qn(t) = fCW > = E fCW " ( 3? ) 

Proof. We check that the right-hand sides satisfy the recurrences in (5) (case a = 1 or 
a = 2), but we only detail the case of R n (case a = 2), the other ones being similar. 
We have R$ = 1, and this corresponds to the snake (0), 1, (0). Also R\ = 2t, and this 
corresponds to (0), 1, —2, (0) and (0), 2, —1, (0). Suppose the result is proved for R n -i, 
then to prove it for R n we distinguish three kinds of elements in (the convention 
is 7To = 7T n +2 = 0). We will denote by tt~ the number 7Tj — 1 if 7Tj > and 7Tj + 1 if 

TTi < 0. 

• First, suppose that |vr n+ i| = 1, hence ir n+ i = (— l) n . Let tt' = tt^ ...tt~. De- 
pending on the parity of n, we have either TT n < ir n+ i = 1 or 7r n > 7r rt+ i = — 1, 
and it follows on one hand that (— l) n vr n < 0, on the other hand that vr n 7r rt+ i < 0. 
Thus tt' G and sc(7r') = sc(7r) — 1. The map tt \— > tt' is bijective, and with the 
recurrence hypothesis, it comes that the set of tt E x with |vr n+ i| = 1 is counted 
by tRn-i- 

• Secondly, suppose that tt\ = 1, and let tt' = —tt^, ■ ■ ■ , — vr~ +1 . From tt\ > TT2, we 
obtain TT2 < 0, hence —tt^ > 0. It follows that tt' £ 5^°. Moreover, since tt\tt2 < 
we obtain sc(7r') = sc(7r) — 1. With the recurrence hypothesis, it comes that the 
set of tt E 5°+! with tt\ = 1 is also counted by tR n -i. 

• Thirdly, suppose that there is j E {2. . . n} such that | vr^- 1 = 1. We have either 
TTj-i > ttj < TTj+i or TTj-i < tt j > TTj+i, and it follows that TTj-i and TTj+i have 
the same sign. We will obtain the term R' n _i in the subcase where ttj have also 
the same sign as TTj-i and vrj+i, and the term t 2 R' n _ 1 otherwise. Let us prove the 
first subcase, the second will follow since it suffices to consider snakes of the first 
subcase where iTj is replaced with — ttj. 

Let 7r' = 7rf > . . . , 7rJ_ 1; —ttJ +1 , . . . , —tt~ +1 . Then we can check that tt h-> {tt' , j) 
is a bijection between tt of the first subcase, and couples (tt 1 , j) where tt' E and 
j E {2 . . . n} is such that tt'a_iTt'a < 0. At the level of generating function, choosing 
a sign change in tt' is done by differentiation with respect to t, and this explains 
why tt of the first subcase are counted by R' n _i- 

Adding the above terms, we obtain (1 + t 2 )R' n _ 1 + 2ti? n _i = R n . This completes the 
recurrence. The difference in the case of Q n (respectively, P n ) is roughly that only the 
second and third points (respectively, only the third point) are to be considered. I 

The exponential generating function of {Qn(t)}n>o can be directly obtained from 
snakes in the form 

oo n 

E„ , s z sec z 

n=0 

To prove this, let tt E 5^, and consider its unique factorization tt = fi ■ ■ ■ fk, where 
each /fc is a maximal non-empty factor such that all its entries have the same sign. 
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Each of these factor has odd length, except possibly the last one. So it is convenient 
to think that if the last factor has odd size, there is also an empty factor fk+i at the 
end, and consider the new factorization tt = f\ . . . /,- (with j = k or j = k + 1). Note 
that sc(7r) = j — 1. The factorization proves (38), because tt is built by assembling 
alternating permutations of odd size /i . . . fj-i, and an alternating permutation of 
even size fj. 

Using the snakes, we can also directly obtain that 



oo 

V 

n=0 \n=0 



7 I 7 \ 



Indeed, let tt = (0), tt\ . . . 7r ra+ i, (0) be among snakes counted by R n (t), i.e. tt £ 5^%- 
There is j G {1 . . . n+1} such that \iTj\ = n+1. We define tt' = (0), tt\ . . . TTj-i, (fTj) and 
tt" = (0), (— l) n 7r n+ i, . . . , (— l) n TTj + i, ((— \) n TTj). After some relabelling, these define 
two snakes in and 5° We have sc(7r) = sc(7r') + sc(7r"), and tt is built by 

assembling tt' and tt" , so this proves (39). 

Another interesting question concerning snakes is the following. For a given per- 
mutation a = o\ . . . a n S 6 n , the problem is to choose signs e = t\ . . . e n £ 
such that tt = ei<7i, . . . , e n a n is a snake. Arnol'd [1] described the possible choices 
such that tt G 5^ in terms of ascent and descent sequences in the permutation. Let 
i G {2 . . . n — 1}, then a case-by-case argument from [1] show that: 

• if Gi-x < Oi < a i+ i, then e« / e i+ i, 

• if Oi-i > Oi > a i+ i, then ej_i ^ e i; (40) 

• if cTj_i > Oi < o i+ i, then = e i+ x. 

Following this, we can answer the problem of building a snake from a permutation in 
terms of some permutations statistics. We can apply this argument to the three sets 
S n , 5°, and 5°°, essentially by varying the conventions on cr(0) and oin + 1). 

Definition 3.5. Let o = o% . . . o n G 6 n , we make the convention that oq = o n +\ = 
n + We also need to consider other conventions, so let 6° and be "copies" of 
the set & n , where: 

• If o G 6°, we set oq = and o n+ \ = n+1. 

• If o G 6°°, we set oq = o n+ i = 0. 

Let o be in one of the sets & n , (3^, or <5^°. An integer i G [n] is a valley of o if 
<Tj_i > Oi < Oi + \, a peak if <Tj_i < Oi > crj+i, a double descent if crj_i > Oi > Oi + \, 
and a double ascent if <7j_i < o"j < o"j + i. Let va(cr) denote the number of valleys in 
o, pk(o") denote the number of peaks in <r, and dda(<r) denote the number of double 
descents and double ascents in o. 

For example, we see 1, 3, 2 G 63 as (0), 1, 3, 2, (4). Then 3 is a valley. On the other 
hand, we see 1, 3, 2 G ©3° as (0), 1, 3, 2, (0). Then 3 is a double descent. 

Proposition 3.6. When n > 1, we have: 

P n (t)= ^ t dda W(i + t 2 ) va W = ^ t dda W(l + t 2 )P k W, (41) 



Q n (t) = Y t dda W(l + t 2 ) va W, R n {t)= t dd < a \l + t 2 )™( a \ (42) 

The last two identities are also true when n = 0. 
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Proof. This follows from Theorem 3.4 and (40), which now can be extended to the 
cases i = 1 or i = n using the conventions on o~q and (7 n +i (note that we also need 
to consider eo and e n +i). Let us first prove the case of P n , a £ S n and ir £ 6 n . 
We have eo = —1 since we require that ttq = — (n + 1). The first two rules in (40) 
show that to each double ascent or double descent we can associate a sign change in 
7r, and that the only remaining choices to be done concern the valleys. At a valley 
(7j_i > Oi < o-j + i, we can have either ej_i = e» = e«+i, or ej_i = — ej = ej+i, i.e. or 2 
sign changes. This proves the first equality in (41), and the second follows by taking 
the complement permutation. The cases of Q n and R n are proved similarly, except 
that we don't consider eo and we have e% = 1. I 

An important property of the statistics va(cr) and dda(cr) is that they can be followed 
through the bijection of Frangon and Viennot (see for example the book [10]). Omitting 
details, this gives combinatorial proof of: 

oo 

Qn{t)z n =2{(2h + l)t, (1 + t 2 )h 2 ) , (43) 

n=0 

oo 

Rn(t)z n = 3( (2/i + 2)t, (1 + t 2 )h(h + 1) ) , (44) 

n=0 

which is the particular case q = 1 of Proposition 2.5. Unfortunately, it seems difficult 
to follow the parameter q through this bijection, and obtain nice statistics on 6° and 
corresponding to Q n (t,q) and R n (t,q). 

3.2. Cycle-alternating permutations. From the combinatorial interpretation in 
terms of snakes, we can derive other ones using simple bijections. We need some 
general definitions concerning signed permutations. Let 7r G S^. 

Definition 3.7. Let neg(7r) be the number of i £ [n] such that 7Tj < 0. For any orbit O 
of 7r, let — O = {—x | x G O}. A cycle of ir is an unordered pair of orbits {Oi, O2} such 
that —0\ = 02- It is called a one-orbit cycle if 0\ = O2 and two-orbit cycle otherwise. 
(This is known to be a natural notion of cycle since there is a decomposition of any 
signed permutation into a product of disjoint cycles.) For any signed permutation ir, 
its arch diagram is defined as follows: draw on the horizontal axis 2n nodes labelled by 
the integers in [rt] in increasing order, then for any i draw an arch from i to such 
that the arch is above the horizontal axis if i < ir(i) and below the axis otherwise. 

See Figure 1 for an example of an arch diagram. Since the labels from — n to n 
are in increasing order we do not need to explicitly write them. See also Figure 2 
further for another example. We did not specify what happens with the fixed points 
in the arch diagram. Although this will not be really important in the sequel, we can 
choose to put a loop ® at each positive fixed point and a loop 5 a t each negative 
fixed point. This way, the arch diagram of a signed permutation is always centrally 
symmetric. 




Figure 1. The arch diagram of tt = 2, -4, 6, 1, -3, 7, 5 
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Definition 3.8. Let tt G ©„. An integer i G [nj is a cycle peak of tt if 7r _1 (i) < i > 
7r(z), and it is a cyc/e valley of 7r if vr _1 (i) > i < n(i). The signed permutation 7r is 
cycle- alternating if every i G [n] is either a cycle peak or a cycle valley. Let C n C 6^ 
denote the subset of cycle-alternating signed permutations, and C° C C n be the subset 
of it G C n with only one cycle. 

Lemma 3.9. Let tt G C°. Then n is even (respectively, odd) if and only if tt has a 
two-orbit cycle (repectively, a one-orbit cycle). 

Proof. In a cycle-alternating permutation, each orbit has even cardinal because there 
is an alternance between cycle peaks and cycle valleys, so if there are two opposite 
orbits of size n, necessarily n is even. 

In the cycle-alternating permutation tt, —n is a cycle valley and n is a cycle peak. 
If there is one orbit, it can be written 7T = (n, ii, . . . , i n -l> — n , ~h, ■ ■ ■ , —in-l) in cycle 
notation. Because of the alternance of cycles peaks and cycle valleys, it follows that 
n is odd. I 

Theorem 3.10. We have : 

Pn{t)= Y, Q^)=E tnegW < Rn ^= E tneg(fT) ' ( 45 ) 

neC° n+1 7reC„ neC° n+2 

7Tl>0 

Proof. There are simple bijections between snakes and cycle-alternating permutations 
to prove this from Theorem 3.4. 

Let us begin with the case of P n , so let tt G S n . Suppose first that n is odd, 
hence ttq = ir n+ i = —(n + 1). We can read tt in cycle notation and consider the 
two-orbit cycle tt' = (ttq, tti, . . . , vr n )(— 7To, — tti, . . . , — n n ). Then %' G C° +1 . When n 
is even, we have — ttq = 7r n+ i = n + 1, and the bijection is defined by taking tt' = 
(ttq,tti, . . . , 7r n , — 7To, — tt\, . . . , — 7r n ). The map tt i — y 7r' is clearly invertible, and using 
the previous lemma, we have a bijection between S n and C° +1 . Let i E {0, . . . , n}, then 
7Tj7Tj + i < is equivalent to 7r r ( | tt^ | ) = — |7Tj + i| < 0, and it follows that sc(7r) = neg(7r'). 
This proves the first equality of the theorem. 

The case of R n is slightly different. In this proof, we will denote by irf the number 
7r.j + 1 if TTi > 0, and Tti — 1 if tt? < 0. Let tt G 5^+1- If n + 1 is odd, hence 7r„ + i > 0, we 
define tt' = (1,^, . . . , 7r^ +1 )(— 1, — irf , . . . , —tt^ +1 ). If n + 1 is even, hence tt„,+i < 0, 
we define tt' = (1, Trf , . . . , vr+ +1 , —1, — tt^ , . . . , — tt^ +1 ). In each case, it is easily checked 
that this defines a bijection from to the subset of tt £ C° +2 satisfying tt\ > 0, so 
that a change of sign in tt correspond to a i > with vr'(z) < 0. 

The case of Q n requires more attention. Let tt S 5^, we need to split it into blocks 
that will correspond to the cycles of an element tt' G C n . This can be done by a 
variant of Foata's fundamental transform [17]. Let a\ < ■ ■ ■ < at be the left-to-right 
maxima of |7q| . . . |7r n | G G n (recall that u is a left-to-right maximum of a G 6 n 
when o(i) < a(u) for any i < u). We consider the factors of tt, f\ = TT ai . . .7r a2 _i, 
/2 = K a2 . . . 7T a3 _i, up to ft = TTa k • • • 7r n (we can take the convention ak+i = n + 1). We 
form tt' by putting together cycles, so that if fi has odd length it gives a one-orbit cycle 
(■7T ai . . . TT ai+1 -i, —TT ai ■ ■ ■ —TTa t+1 -i), and if it has even length it gives a two-orbit cycle 
(■K at . . . TT ai+1 -i)(—TT ai . . . — 7r ai+1 _i). The definition of aj shows that n ai is greater than 
7r ai+ i . . . 7r ai+1 _i, and consequently these cycles are indeed alternating. The inverse 
bijection is easily deduced: we can write tt' G C n as a product of cycles tt' = c\ . . . cp., 
where the cycles are sorted so that their maximal entries are increasing. To each q we 
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associate a word fi so that Cj is either (/j, — /j) or (fi)(—fi), and such that the first letter 
of fi is the maximum entry of Cj. Then we can form the snake tt = /i, €2/2, ■ ■ ■ , ejfc/fc, 
where the signs 62, ... , £fc G {±1} are the unique ones such that tt is alternating. 

It is in order to give examples of the previous two bijections. If we start from 
a snake tt = (0), 4, —2, —1, —5, 3, —6, (7) G ^6, the left-to-right maxima of \ir\ are 
1,4,6, and this gives tt' = (4, -2, -1, -4, 2, l)(-5, 3)(5, -3)(6, -6) G C 6 . As for the 
other direction, let vr' = (5, -2, 3, -5, 2, -3)(6, -l)(-6, 1)(7, -4)(-7, 4) G C 7 . It is a 
product of three cycles so that the snake tt is formed by putting together three words 
5, —2, 3, and 6, —1, and 7, —4. We need to change the signs of the second and third 
words to obtain a snake, which is tt = (0), 5, —2, 3, —6, 1, —7, 4, (—8). 

One can check on the examples that in both cases we have sc(7r) = neg(7r'), but 
this does not immediately follows the construction and need to be proved now. Let 
a, b > 0, such that vr'(a) = —b, we can associate a sign change in tt to each such 
pair (a, b). The two integers are in the same cycle q. If the word fi is of the form 
. . . , a, —6, ... or ... , —a, b, . . . , then either the factor a, —b or —a, b appears in the 
word eifi, and consequently appears in tt as a sign change. The other possibility is 
that fi = b,..., —a which only occur when fi has even length, and fi = b, . . . ,a which 
only occur when fi has odd length. If b (respectively —b) appear in tt, it is greater 
(respectively, smaller) than his neighbor entries, and after examining a few cases it 
comes that the same is true for a. It follows that there is a sign change in tt between 
the last entry of eifi and its right neighbor. This completes the proof. I 

Using cycle-alternating permutations, there is a bijective proof of the fact that 
P n +i = (1 + t 2 )R n . This can be written 

^ £ neg W = (l + i 2 ) ^ t ncg(<T) - (46) 

0-1 >o 

Now, consider the conjugation by the transposition (1, —1), i.e. the map that sends tt 
to (1, — l)vr(l, —1). This is a fixed-point-free involution on C° +2 which proves combi- 
natorially (46), as is easily seen on the arch diagrams. Indeed, in this representation 
the involution exchanges the dot labelled 1 with the dot labelled —1, see Figure 2. We 
have 7r(l) > if and only if its image tt' = (1, — l)vr(l, —1) does not satisfy tt'(1) > 0, 
and if tt(1) > there holds neg(7r') = neg(7r) + 2. This proves (46). Note that it also 
is possible to prove P n +i = (1 + t 2 )R n on the snakes, using Proposition 3.6, but this 
is more tedious. 



Figure 2. The arch diagram of the cycle-alternating permutation 
tt = (1, 4, —5, —2, —3, —1, —4, 5, 2, 3) G C5, and of its conjugate 
(l,-l)7r(l,-l). 

Using the cycle-alternating permutations, we can also obtain the exponential gen- 
erating function of {P n }n>o, knowing the one of {Q n }n>o- Indeed, a cycle-alternating 
permutation is an assembly of cycle-alternating cycles, in the sense of combinatorial 



ENUMERATION OF SNAKES AND CYCLE-ALTERNATING PERMUTATIONS 



13 



species, and this gives immediately: 

y n+l 



z 



(n + 1)! 



n=0 

and after differentiation: 



\n=0 



t 



cos z — t sin z 



ra=0 



??,! 



sin z + t cos z 
cos z — t sin z 



(47) 



(48) 



3.3. The g-analogs. It is possible to give combinatorial interpretations of Q n (t, q) 
and R n (t, q) using cycle-alternating permutations, and the notion of crossing for signed 
permutations which was defined in [5]. But in the case of cycle-alternating permuta- 
tions, there is a more simple equivalent definition which is the following. 

Definition 3.11. A crossing of ir £ C n is a pair (i,j) E [nj 2 such that i < j < ir(i) < 
7r(j). Equivalently, it is the number of intersections of two arches above the horizontal 
axis in the arch diagram of ir. We denote by cr(-7r) the numbers of crossings in ir. 

Because of the symmetry in the arch diagram, cr(7r) is also half the total number 
of intersections between two arches. For example, the two permutations in Figure 2 
have respectively 2 and 3 crossings. This kind of statistic can be followed through 
bijections between permutations and paths [4], but let us first show that Q n (t,q) and 
R n (t,q) are related with some weighted Dyck prefixes (a Dyck prefix being similar to 
a Dyck path except that the final height can be non-zero). 

Definition 3.12. Let V n be the set of weighted Dyck prefixes of length n, where 

• each step between heights h and h+1 has a weight q l for some i € {0, . . . h+1}, 

• each step \ between heights h and h + 1 has a weight q % for some i 6 {0, . . . h}. 

Let V' n C V n be the subset of paths p such that there is no step /* starting at height 
h with the maximal weight q h+l . For any p G V n , let fh(p) be its final height, and let 
w(p) be its total weight, i.e. the product of the weights of each step. 



Proposition 3.13. We have: 

Qn(t,q)= £ t^q w(p) , 
P eV' n 



Rn(t,q) 



w(p) 



(49) 



Proof. In the proof of Proposition 2.5, we have used the matrices Z?2 + U2D2U2 and 
D2 + D2U2U2, thought of as transfer-matrices for walks in the non-negative integers. 
We can do something similar with the matrices 



D 1 + U x DxUx 



(0 


[1] 9 


(0)\ 


[1], 





[2], 






'•. 


\(0) 




'••/ 



Dx + DxUxUi 



(0 


[1] 9 


(0)\ 


[2], 


[2], 






[3], 




\(0) 




•7 



(50) 



Recall that we have (20), where W\ = (1, t, t 2 , . . . ). In terms of paths, the coefficients 
in the matrices and vectors have the following meaning: from V\ = (1, 0, . . . )*, we only 
consider paths starting at height 0, and from W\ = (l,t,t , . . . ) there is the weight 
t^ip) for each path p. The coefficients (h + l,h) in the matrices give the possible 
weights on steps from height h to h + 1, the coefficients (h, h + 1) in the matrices 
give the possible weights on steps \ from height h + 1 to h. The result follows. I 
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Proposition 3.14. There is a bijection ^ : C n 
cr(7r) = w(p) if the image of tt is p. 



V' n such that neg(7r) = fh(p) and 



Proof. We can use the the bijection ^ fz from [4] between permutations and weighted 
Motzkin paths. To do this, we identify to a subset of &2n via the order-preserving 
bijection \n\ — > [2n]. This subset is characterized by the fact that the arrow diagrams 
are centrally symmetric, and via ^fz from [4], it is in bijection with some weighted 
Motzkin paths that are vertically symmetric. Cycle-alternating permutations corre- 
pond to the case where there is no horizontal step, i.e. they are in bijection with some 
weighted Dyck paths that are vertically symmetric. Keeping the first half of these 
Dyck paths gives the desired bijection. I 

For convenience, we rephrase here explicitly the bijection of the previous proposition. 
Let 7r G C n , then we define a ^(ir) G V' n such that, for any j G {— n . . . —1}: 

• the (n + 1 + j)th step is /* if vr _1 (j) > j < ir(j) and \ if vr _1 (j) < j > vr(j'), 

• if the (n + 1 + j)th step is /*, it has a weight q k where k is the number of i such 
that is a crossing, i.e. i < j < 7Tj < Wj, 

• if the (n + 1 + j)th step is \, it has a weight q k where k is the number of i such 
that (— i, — j) is a crossing, i.e. —i < —j < — 7Tj < — 7Tj. 

See Figure 3 for an example. 












to 

















































Figure 3. The bijection ^ from C n to T 3 ^ i n the case of tt 
3,4,-7,2,6,1,5. 



Through this bijection we obtain combinatorial models of Q n (t, q) and R n (t, q) using 
the notion of crossing. 

Proposition 3.15. We have 

Q n (t,q)= i? n (t,g)= ^ t n egW -yrM (gl) 

1"n + l<0 

Alternatively, the condition vr n+ i < in the second identity can be replaced with 
7r" 1 (n+ 1) < 0. 

Proof. The first identity follows from Proposition 3.13 and the bijection ^ between C n 
and V' n , but the second one is not as immediate. 

First, note that the bijection tt i— > tt" 1 stabilizes the set C n and preserves the number 
of crossings, since it is just an horizontal symmetry on the arch diagrams. This proves 
the fact that we can replace ir n+ i < with 7r _1 (n + 1) < 0. 

Then, let us consider the subset V^ +1 C V' n+1 of elements p such that there is no 
step \ from height h + 1 to h with weight q h (note that this condition implies there 
is no return to height 0). There is an obvious bijection from to V n , because 

removing the first step in p G "P^'+i §i yes a P a th which is in V n with respect to the 
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shifted origin (1,1). This proves: 



Rn(t,q)= tHP) ~ 1( l 



, w (p) 



(52) 



It remains to show that the image of V^+x y i a ^ is precisely the set of it G C n +i 
such that 7r _1 (n + 1) < 0. Essentially, this follows from the fact that there are some 
steps in the path which characterize the right-to-left minima in it, and we can use 
Lemma 3.2.2 from [14]. In our case, the subset V" +1 C V' n+1 corresponds to a subset 
of &2n+2 (via the identification 6^ +1 — > &2n+2), more precisely to the subset of a 
having no right-to-left minima among 1 . . . n + 1, i.e. <7 -1 (l) > n + 1. In terms of 
7r G &n+i> this precisely means that ir~ 1 (—n — 1) > and completes the proof. I 

Remark 3.16. Independently from our work, Chen, Fan and Jia [3] gave a bijection 
between snakes of type B and Dyck prefixes such that there are h + 1 possible choices 
on each step between height h and h + 1, and their statistic a(n) (see Theorem 4.5) 
is the same as our statistic sc(-7r), though defined differently. We can have such a 
bijection by composing two we have presented, but their bijection is not the same as 
ours and does not seem to be directly related. To advocate our bijection S n — > V' n , 
note that going through the set C n we just had to adapt two known bijections: Foata's 
fundamental transform [17] (in the direction S n — > C n ), then Corteel's bijection ^ fz 
[4] originally given by Foata and Zeilberger (in the direction C n — > V' n ). 

3.4. Another variant of ^fz- With another adaptation of this bijection, we can 
link C n and some weighted Motzkin path, so that we obtain a combinatorial proof of 
a J- fraction generalizing (43) and (44) in Theorem 3.17 below. As a consequence of 
Theorem 3.10, we have: 



where cyc(7r) is the number of cycles, indeed since there is a "cycle structure" in 
C n it suffices to take the ath power of the exponential generating function to have 
a parameter a counting cyc(-7r). With the recurrence for Q^h\t) in (5) and with the 
same method as in the case of Q n and R n , we can prove a continued fraction for the 
ordinary generating function of {Q^ (t)} n >o- We prove here bijectively this continued 
fraction, by going from C n to weighted Motzkin paths. To this end, we use another 
representation of a signed permutation, the signed arch diagram. Let it G S^, this 
diagram is obtained the following way: draw n dots labelled by 1 . . . n from left to 
right on the horizontal axis, then for each i G [n], draw an arch from i to | vr^ | , and 
label this arrow with a + if 7Tj > and with a — if tti < 0. We understand that as in 
the previous arch notation, the arch is above the axis if i < | vr^ | and below otherwise. 
See the left part of Figure 4 for an example. (In the case | vr^ | = i, the arch is a loop 
attached to the dot i.) A case-by-case study shows that tt is cycle-alternating if and 
only if it avoids certain configurations in the signed arch diagram, which are listed in 
the right part of Figure 4 (± being either + or — ). These forbidden configurations are 
of course reminiscent of (40). 

Once we know the forbidden configurations, it is possible to encode the signed arch 
diagram of it G C n by a weighted Motzkin path of n steps, in a way similar to ^fz 
from [4]. The diagram of tt is scanned from left to right, and the path is built such 
that: 




(53) 
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+ 



Figure 4. The signed arch diagram of tt = 2, —4, —6, 5, —3, —1,7, and 
the forbidden configurations in the signed arch diagram of a cycle- 
alternating permutation. 

• If the ith node in 7r is <^ then the ith step of p if /*, moreover this step has label 
+ (resp. — ) if the two arches starting from the ith node have label + (resp. — ). 

• If the ith node in ir is ^then the ith step of p if — K Moreover this step has label 
(j) if the left strand in the ith node is connected to the jth strand of the "partial" 
signed arch diagram (say, from bottom to top). See for example the left part of 
Figure 5 where there are three possible choices to connect this kind of node. We 
can see that the possible labels are (1), . . . , (h) where h is the starting height of 
the ith step in the path. 

• Similarly if the ith node in tt is ^then the ith step of p if — >, and this step has 
label (j') if the ith node is connected on the left to the jth. strand of the "partial" 
signed arch diagram (say, from top to bottom) . 

• If the ith node in tt is 8 then the ith step is — > with label (0). 

• If the ith node in tt is )• then the ith step is \ with a label (j, k) where j and 
k respectively encode where the upper and lower strands are connected, as in the 

case of ^and ^\ above. 



+ 



+ 



















\/ 




+ 


(1) 


/ 






/ 











+ 



Figure 5. Construction of the bijection from the signed arch diagram 
of 7r G C n to a labelled Motzkin path. 



Quite a few bijections of this kind are known so we will not give further details. As 
for the parameters a and t, they are taken into account the following way: 

• For each arch with label — linking i and j where i < j, we give a weight t to the 
ith step in the path. Thus there is a weight t 2 on the step /* with label — , and a 
weight t on each step — >. 

• For each i which is the maximal element of a cycle, we give a weight a to the ith 
step. Thus each step —> with weight (0) has a weight a, since these correspond to 

nodes 8 or cycles of the form (i, — i). Besides, there is a weight a on some of the 
step \. More precisely, there are h 2 possible labels on a step \ starting at height 
h, and h of them have this weight a. Indeed, for each choice of where to connect 
the upper strand of the node ^> there is a unique choice of where to connect the 
lower strand so that a new cycle is "closed". 
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All in all, using the classical correspondence between weighted Motzkin paths and 
J-fractions, this bijection proves the following theorem. 

Theorem 3.17. Let Qn\t) be the polynomials in a and t with exponential generating 
function (cos z— t sin z)~ a , or defined combinatorially by (53). Then, using the notation 
in Definition 2.4, we have: 

oo 

Y^Q^(t)z n = Z{ (2h + a)t, h(h-l + a)(l+t 2 ) ). (54) 

n=0 

Apart the combinatorial side, this result was already known by Stieltjes [22]. Besides 
his analytical proof, it is also possible to obtain the continued fraction by showing that 
the exponential generating function satisfies some particular kind of addition formula 
which permits to make use of a theorem of Stieltjes and Rogers, see for example 
[10, Theorem 5.2.10]. We thank Philippe Flajolet for communicating these facts and 
references. 



4. Increasing trees and forests 

So far, we have only used recurrence relations and bijections to derive our new 
models of the derivative polynomials. Some more elaborate methods apply pretty well 
in our case to give combinatorial models in terms of increasing trees and forests. We 
first present how to use the combinatorial theory of differential equations, as exposed 
by Leroux and Viennot in [16], and secondly how to use recent results of Blasiak and 
Flajolet related with operators and normal ordering [8]. 

4.1. Increasing trees via differential equations. An archetype example in the 
combinatorial theory of differential equation is the one of the tangent and secant 
functions (see [16]), and it has become a classical method to show that Euler numbers 
E n count some increasing trees, i.e. labelled trees where labels are increasing from the 
root to the leaves. Here we have a system of two differential equations similar to the 
one of tangent and secant, except that an initial condition is given by the parameter t 
instead of 0. This will give rise to more general trees where we allow some leaves with 
no label having a weight t. 

Lemma 4.1. Let f = ^2 P n z n /n\ and g = Y2Qn zU / n ^> then we have: 

Jf = l + / 2 /(0)=t, 

\g> = fg g(0) = 1. (55) 

Proof. Of course, this can be checked on the closed form given in (7). Also, /' = 1 + / 2 
can be directly checked on the snakes as in the case of alternating permutations [21], 
and g' = fg follows from / = (log g)' previously seen. I 

It is adequate to rewrite the equations in the following way: 

l f= l + + n Sf2 ' < 56 > 

[9 = 1 + J fd- 
Let us begin with the case of f(z). From / = t + z + j f 2 , and proceeding as in 
[16], f(z) counts increasing trees that are recursively produced by the following rules, 
starting from an isolated node marked by /: 
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• a node marked by / can become a leaf with no label (this corresponds to the term 
t in / = t + z + J*/ 2 , f(z), so these leaves will have a weight t), 

• a node marked by / can become a leaf with an integer label, this label being the 
smallest integer that does not already appear in the tree (this corresponds to the 
term z in f = t + z + f f 2 ), 

• a node marked by / can become an internal node having an integer label and 
two (ordered) sons marked by / (this corresponds to the last term f f 2 in / = 
t + z + f f 2 ). As before, the integer label is the smallest integer that does not 
already appear in the tree. 

More precisely, the coefficient of z n in f(z) counts these trees having n integer labels. 
See Figure 6 for an example of tree produced by these rules. 



© 




Figure 6. Tree produced via the equation / = t + z + J f 2 . 

The case of g(z) is quite similar. Starting from an isolated node marked by g, each 
node marked by g can become either a leaf with no label, or an internal node having 
two sons respectively marked by / and g. Note that we need the production rules for 
/ to build a tree counted by g, and note also that g can produce an empty leaf which 
has no weight t, contrary to the empty leaves produced by /. 

The two kinds of tree can also be given a non-recursive definition. 

Definition 4.2. Let T n be the set of complete binary trees, such that: 

• except some leaves that are empty, the nodes are labelled by integers so that each 
i £ [n] appears exactly once, 

• labels are increasing from the root to the leaves. 

Let em(T) be the number of empty leaves of T G T n , and let T* C T n be the subset 
of trees such that the rightmost leaf is empty. 

The production rules for / (respectively, g) can be checked on the set T n (respectively 
T*), so that the result of the above discussion is the following. 

Theorem 4.3. We have: 

Pni t) = y, * em(T) , Quit) = * em( ™- ( 57 ) 

TeT n TeT* 

4.2. Increasing trees via normal ordering. Some recent results of Blasiak and 
Flajolet [2] directly apply to this context and also give models of P n (t) and Qn (t) in 
terms of increasing trees and forests. Let D be the derivation with respect to t and 
U the multiplication by t (i.e. we are in the particular q = 1 of operators defined 
in (11)). There holds DU — UD = I. A general idea in [2] is that the coefficients 
Cij in the normal form of f(D,U) as in (17), at least for certain particular forms 
of f(D,U), naturally counts some labelled directed graphs which are produced by 
connecting some "gates". In the present case, when f(D,U) is (D + UDU) n for P n 
and (D + all + UUD) n for and we obtain some increasing trees and forests. 
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The theorem of this subsection is a direct application of a main result of Flajolet and 
Blasiak [2, Theorem 1], so as a proof we will roughly explain some ideas leading to the 
definition below and refer to [2] for more details. 

Definition 4.4. Let T n be the set of plane rooted forests satisfying the following 
conditions: 

• each root has exactly one child, and each of the other internal nodes has exactly 
two (ordered) children, 

• there are n nodes labelled by integers from 1 to n, but some leaves can be non- 
labelled (these are called empty leaves), and labels are increasing from each root 
down to the leaves. 

Note that the trees forming a forest are unordered. Let lA n C T n be the subset of 
trees, i.e. forests with one connected component. For any tree or forest T, let em(T) 
be the number of empty leaves, and let cc(T) be its number of connected components. 



For example, there are 11 elements in T-& and they are: 




Theorem 4.5. We have: 

P n (t)= * em(T) , Q { n\t)= Y, a cc(T) i em ( T ). (58) 

For example, the forests in .F3 given above illustrate Q3 = 6i 3 + 5t. The last four 
elements of the list are the trees, and illustrate P2 = 2i 3 + It. By counting with a 
weight 2 on each connected component, we obtain from this list that R3 = 24i 3 + 16t. 

Proof. From the definition of derivative polynomials in terms of D and U, and the 
relation DU — UD = I, we have: 

P n (t) = {D + UUD) n t, Q { n\t) = (D + aU + UUD) n l. (59) 

Let us consider the case of Q { n ] . Let f n (D, U) = (D + aU + UUD) n and CrMj - the 
coefficient of U l D 3 in its normal form, as in (17). Then from [2, Theorem 1], c n ^^j 
counts labelled diagrams obtained by connecting three kinds of "gates", one for each 
term in D + all + UUD. More precisely to each term U k D e we associate a gate 
consisting of one node with k outgoing strands and £ ingoing strands. See the left part 
of Figure 7. Then, Cn^j count labelled diagrams obtained by connecting n of these 
gates such that: 

• i outgoing strands and j ingoing strands are not connected, 

• all other strands are connected, so that each ingoing strand is connected with an 
outgoing strand and these form a directed edge, 

• the gates labelled by the integers in {1 . . . n}, and labels are increasing when we 
follow a directed edge, 
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• at each node, the ingoing strands on one side and the outgoing strands on another 
side are ordered, 

• there is a weight a at each gate corresponding to the term U. 

We have 

QW(t) = (D + aU + UUD) n l = ( a,j uiDj ) 1 = Ci > of ( 6 °) 

i,j>0 i>0 

so that we can only consider diagrams with no unconnected ingoing strand, and t 
counts the unconnected outgoing strands. 

The labelled diagrams described by the above rules are essentially the same as 
elements in T n : it suffices to add an empty leaf at each unconnected outgoing strand 
to see the equivalence. It is clear that the node corresponding to the term U will 
appear exactly once in each connected component of the labelled diagrams, so the 
parameter a counts indeed the connected components. 

In the case of P n , we can also consider f n (D,U) = (D + UUD) n U and c n< ij the 
coefficient of U l D 3 in its normal form, as in (17). This case is somewhat different 
since f n (D, U) is not the nth power of some expression, but similar arguments apply 
as well: the labelled diagrams that appear have n + 1 gates, the gate labelled 1 is of 
type U, all other gates are of type D or UUD. These labelled diagrams are the same 
as elements in U n +i, as in the previous case we just have to add an empty leaf to each 
unconnected outgoing strand to see the equivalence. As previously said, we refer to 
[2] for precisions about this proof. I 




Figure 7. The "gates" corresponding to terms in D + all + DUU, 
and a labelled diagram obtained by connecting them. 

There is a simple bijection between T n and U n +i' given T S T n , relabel the nodes by 
i i— > i + 1, then add a new node with label 1 on top of the root. There is also a simple 
bijection between T* and T n : let T £ T*, remove the rightmost leaf, as well as all 
edges in the path from the root to the rightmost leaf, then the remaining components 
form the desired forest. See Figure 8 for an example. 




Figure 8. The bijection from 7^* to J 7 , 



Thus, Theorems 4.3 and 4.5 are essentially equivalent although obtained by different 
methods. It is also in order to give a bijection between U n and C°, and by applying this 
bijection componentwise it will give a bijection between J- n and C n . It is more practical 
to give the bijection from Tn to S n (recall that we already have simple bijections 
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T n — > U n +i and S n — > C° +1 ), so let T £ T n - Consider the "reading word" of this tree 
(it can be defined by w(T) = w(Ti) i w{T2) if the tree T has a root labelled i, left son 
T\ and right son T2). This word contains integers from 1 to n, and some letters (say 
o) to indicate the empty leaves. The first step is to replace each i with n + 1 — i in this 
word. To obtain the snake, replace each integer i by (— 1) J+ i where j is the number 
of o before i in the word, then remove all o. See Figure 9 for an example. 




3,2,o,l,o,4,o -> 2,3,o,4,o,l,o -> (-5), -2, -3, 4, -1,(5) 



Figure 9. The bijection from T n to S n . 

Actually, this bijection can be highlighted by Proposition 3.6 and the discussion 
leading to it. Indeed, it is just a variant of a classical bijection between permutations 
and unary-binary increasing trees [23] , where double ascents and double descents cor- 
respond to nodes having only one child, and valleys correspond to leaves. Via these 
bijections, removing the empty leaves of T £ T n is the same as taking the absolute 
value of a snake S £ S n . In the other direction, T seen as a unary-binary tree together 
with the data of the empty leaves, is the same as a permutation together with a choice 
of signs making it into a snake. 

Conclusion 

There is still a wide range of results on alternating permutations that might have a 
counterpart in the case of snakes, for example one could ask if Andre permutations and 
simsun permutations [21], both counted by E n , have natural generalizations in signed 
permutations (as a first answer, see [18] for a definition of Andre signed permutation, 
counted by S n ). Another problem is to give more combinatorial meaning to the Arnold- 
Seidel triangles [1, 6] used to compute the integers S n (see [11] for recent work on this 
subject in the case of integers E n ). As for g-analogs, we have generalized the numbers 
E n (q) defined in [12] but there is another q-analog of E n related with inversions and 
major index [21] that might also be extended to snakes. 

Though we only discussed the type B ones, ArnoPd [1] also defined some type D 
snakes counted by the Springer number K(D n ). Besides, Hoffman showed that K{D n ) 
is equal to -P n (l) — QniX)-> so that it is also the number of tt £ S n such that 7ri < 0. As 
usual with type D, it is much more difficult to obtain enumerative or bijective results 
with either of these two families, they are not even a subset of the group of even- 
signed permutation D n as one could expect. Having said that, one possible direction 
towards a better understanding of these objects could be the geometric definition of 
the Springer number K(W) in terms of hyperplane arrangements [1]. For example, 
in this context it would be very interesting to give the polynomial Q n (t) a geometric 
meaning that would refine Q n (l) = K(<S%). 
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